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EFFECTS  OF  SHEAR  DEFORMATION  AND  ANISOTROPY 
ON  THE  THERMAL  BENDING  OF  LAYERED  COMPOSITE  PLATES 


J.  N.  Reddy  and  Y.  S.  Hsu 

School  of  Aerospace,  Mechanical  and  Nuclear  Engineering 
University  of  Oklahoma,  Norman,  OK  73019 

A  finite-element  formulation  of  equations  governing  layered  anisotropic 
composite  plates  subjected  to  thermal  and  mechanical  loadings  is  presented.  An 
exact  closed-form  solution  is  also  presented  for  simply  supported  rectangular 
cross-ply  laminated  plates  under  sinusoidal  loading  to  validate  the  finite 
element  developed  herein.  The  finite-element  results  are  in  good  agreement 
with  the  closed-form  results  and  with  the  results  of  others.  Material  pro¬ 
perties  typical  of  advanced  fiber-reinforced  composites  are  used  to  show 
the  parametric  effects  of  plate  aspect  ratio,  side-to-thickness  ratio,  orien¬ 
tation  of  layers,  and  edge  conditions  on  the  deflections  and  stresses. 

INTRODUCTION 

With  the  increased  use  of  fi ber -reinforced  composites  in  aerospace  and 
mechanical  engineering  structures,  studies  involving  the  thermomechanical 
behavior  of  composite-material  plates  and  shells  are  receiving  greater  atten¬ 
tion.  Most  of  the  previous  research  in  the  field  of  composites  deals  with 
isothermal  problems.  However,  use  of  composites  in  environments  with  large 
temperature  changes  (e.g.,  space  shuttle)  requires  the  knowledge  of  thermally 
induced  deflections  and  stresses.  Further,  thermal  stresses  are  also  induced 
during  the  fabrication  of  composite  materials. 

The  problem  of  thermal  bending  of  anisotropic  plates  was  studied  first 
by  Pell  [1],  who  derived  the  equations  governing  the  transverse  deflection 
of  a  thin  plate.  Generalization  of  Pell's  work  to  heterogeneous  plates 
subjected  to  arbitrary  three-dimensional  temperature  distribution  is  due  to 
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Stavsky  [2],  Recent  studies  in  the  analysis  of  plates  laminated  of  fiber- 
reinforced  materials  indicate  that  the  thickness  effect  (i.e.  shear  deformation) 
on  the  behavior  of  the  plate  is  more  pronounced  than  in  isotropic  plates  [3]. 

The  shear  deformation  theory  that  has  been  proven  to  be  adequate  in  predicting 
the  overall  response  of  laminated  anisotropic  plates  is  due  to  Yang,  Norris 
and  Stavsky  [4],  Based  on  the  Yang-Norris-Stavsky  (YNS)  theory,  Reddy  [5,6] 
developed  a  finite-element  model  that  is  algebraically  simpler  than  previously 
developed  finite  elements  [7-10],  and  yet  possesses  competitive  accuracy. 

The  present  investigation  is  concerned  with  the  application  of  the 
penalty  finite  element  [6]  to  the  thermal  stress  analysis  of  layered  aniso¬ 
tropic  composite  plates.  To  illustrate  the  accuracy  of  the  present  element, 
closed-form  solutions  are  also  presented  for  the  equations  governing  (i.e., 
the  YNS  theory)  simply  supported,  rectangular ,  cross-ply  plates  under  sinusoid¬ 
al  mechanical  and/or  thermal  loadings.  Finite-element  solutions  are  presented 
to  show  the  effects  of  variations  in  geometry,  lamination  parameters,  boundary 
conditions,  and  loading  on  the  shear  deformation  and  thermal  response  of 
statically  loaded  layered  anisotropic  composite  plates. 

GOVERNING  EQUATIONS 

Consider  a  plate  of  constant  thickness t composed  of  a  finite  number 
of  anisotropic  layers  with  arbitrary  orientations.  The  coordinate  system 
is  such  that  the  middle  plane  of  the  plate  coincides  with  the  x-y  plane, 
and  the  z-axis  is  normal  to  the  middle  plane,  1R. 
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where  u,  v,  and  w  are  the  displacements  along  x,  y,  and  z  directions, 
respectively,  uQ  and  vQ  are  the  in-plane  displacements  of  the  mid-plane,  and 
ip  and  ip  are  the  shear  rotations. 

x  y 

The  equilibrium  equations  associated  with  the  YNS  theory  are 


N1,x  *  N6,y  =  pl  •  N6,x  *  NZ,y  :  p2  •  Q1  ,x  *  °2,y  *  '<> 


M,  +  Mc  -  Q,  -  0  ,  M,  +  M„  -  Q„  =  0 
l,x  6,y  yl  6,x  2,y  w2 


(2) 


where  Ni  x  =  al^./ax  etc.,  p^  and  p2  are  the  in-plane  distributed  forces, 
q  is  the  transversely  distributed  load,  and  Ni ,  Qi  and  Mi  are  the  stress 
and  moment  resultants  defined  by 


t/2 

(l,z)o-j  dz  ,  (Q,  ,Q?) 
-t/2 


t/2 


-t/2 


(aXz»Oyz)  dz 


(3) 


Here  a. (i=l  ,2,6)  denote  the  in-plane  stress  components  (a1  =  ax>  a2  =  ay, 


Assuming  monoclinic  behavior  ( i . e . ,  one  plane  of  elastic  symmetry)  for  each 
layer,  the  constitutive  equations  for  an  arbitrarily  laminated  plate  are 
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The  plate  stiffnesses  A^. ,  B^.,  and  are  given  by 
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m+1 


<WDu>-J 


qj^O.z.z2)  dz  ,  (i,j  =  1,2,6) 


m+1 


(5) 


A.  .  =  I 
U  m 
m 


,(m) 


Wlj  dZ  ■  <“  ‘  *-'•  6  '  6‘ji  i>J  '  4'5> 


.(m) 


where  QV.'  are  the  stiffness  coefficients  of  the  m-th  layer  in  the  plate  coordinates, 

*  J 

and  zm  is  the  distance  from  the  mid-plane  to  the  lower  surface  of  the  m-th 
layer.  The  stress  and  moment  resultants,  nI  and  mT,  due  to  thermal  loading 
are  defined  by 


(N{,M I)  =  E  ^  a .(T  ,zT.)  dz  ,  (i,j  =  1,2,6) 

m  Jz  j  J  J 
m 


(6) 


where  are  the  thermal  coefficients  of  expansion  in  the  plate  coordinates, 
and  T  is  the  temperature  change  from  a  reference  state, 

T(x,y,z)  =  TQ(x,y)  +  zT^x.y)  (7) 


Note  that  the  temperature  variation  through  the  thickness  is  assumed  to  be 
linear,  consistent  with  the  plate  theory. 

Substituting  Eq.  (4)  into  Eq.  (2),  we  obtain  the  following  operator 
equation, 

[L]{5}  =  {f}  (8) 

where  {6}  =  {u0,vo>w,4)x,i((y}T,  [L]  is  the  (symmetric)  matrix  of  differential 
operators, 

L11  =  Alldll  +  2A16d12  +  A66d22 

L12  =  (A12+A66^d12  +  A1 6dl 1  +  A26d22  ’  L13  =  0  ’ 

L14  =  Blldll  +  2B16d12  +  B66d22 

L15  =  (B12+B66)d12  +  B1 6dl 1  +  B26d22  =  L24  ’ 

L22  =  2A26d12  +  A22d22  +  A66dll  ’  L23  =  0  » 

L25  =  2B26d12  +  B22d22  +  B66dll 
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L33  =  ■A44dir2A45d12‘A55d22’  L34  =  'A44dl"A45d2 
L35  =  ■A45dTA55d2’  L44  =  °lldll+2D16d12+D66d22'A44 

L45  =  ^D12+D66^d12+D16dn+D26d22-A45’  L55  =  2D26d12+022d22+066dll  ^ 


and  the  components  of  the  generalized  force  vector,  { f } ,  are  given  by 

fl  ■  N(,x  *  N6.y  *  P1  •  f2  ■  "J.x  +  "L  *  p2 

f,  =  q  ,  f .  =  m!  +  Ml  ,  fc  =  m!  +  M. 

3  4  1  ,x  6,y  5  6,x  2,y 


(10) 


In  Eq.  (9),  d.  .  denote  the  differential  operators 

*  J 

=  • «,  -"fo-sfr  -0,1,2) 

I  J 


EXACT  CLOSED-FORM  SOLUTION 


The  boundary-value  problem  associated  with  the  equilibrium  of  layered 

anisotropic  composite  plates  involves  solving  the  operator  equation  (8) 

subjected  to  a  given  set  of  boundary  conditions.  It  is  not  possible  to 

construct  exact  solutions  to  Eq.  (8)  when  the  plate  is  of  arbitrary  geometry, 

constructed  of  arbitrarily-oriented  layers,  and  subjected  to  an  arbitrary  loading 

or  boundary  conditions.  However,  an  exact  closed-form  solution  to  Eq.  (8)  can 

be  constructed  when  the  plate  is  of  rectangular  geometry  with  the  following 

edge  conditions,  loading,  and  plate  construction. 

Boundary  conditions  (freely  supported) 

u(x,0)  =  u(x,b)  =  0,  N2(x,0)  =  N2(x,b)  =0  ^ 

v(0,y)  =  v(a,y)  =  0,  N^O.y)  *  N^a.y)  =  0 

w(x,0)  =  w(x,b)  *  w(0,y)  *  w(a,y)  =  0  (11)  b 


4<x(x,0)  =  ipx( x ,b)  =  0,  M2(x,0)  *  M2(x,b)  =  0 
*y(0,y)  *  tpy(a ,y)  =  0,  M^O.y)  =  M^a.y)  =  0 


f 

liL 
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Loading  (sinusoidal) 

q=q™1  sin  ox  sin  By,  T  “T™1  sin  ax  sin  By,  T1=T,?n  sin  ax  sin  sy 

0  0  0  I  I  ^  i  • 

P-|=P^n  sin  ax  cos  By  ,  P2=P^n  cos  ax  sin  By  ,  a  =  rrnr/a  ,  B  =  mr/b 
and  m  and  n  are  integers. 

Plate  construction  (cross-ply,  i.e.,  em  should  be  either  0°  or  90°) 

A16  =  A26  =  A45  *  0  ’  B16  3  B26  =  0  ’  °16  =  °26  =  0  ’  a6  =  0 

Under  these  specific  conditions,  the  solution  (u.  ,v  ,w,  ^  ^  )  to  Eq.  (8) 

oo  a  y 

is  of  the  form, 

u„  =  Um„  cos  ax  sin  By  ,  v  =  sin  ax  cos  By 
o  mn  J  o  mn 

w  *  W  sin  ax  sin  By  (1< 

\  =  Xmn  cos  ax  sin  »  *y  =  Ymn  sin  ax  cos  &y 

where  Umn,  Vmn>  etc.  are  parameters  to  be  determined  subjected  to  the  con¬ 

dition  that  the  solution  in  Eq.  (14)  satisfies  the  operator  equation  (8). 
Substituting  Eq.  (14)  into  Eq.  (8),  we  get 

[C] {A}  =  {F}  (1‘ 


where 


tA1  =  <v  v  v  vT-  <f>  ■  (iT.fr.Fr.fr.FT)7 


and  the  elements  of  the  coefficient  matrix,  [C],  are  given  by 

C11  =  'Alla  "  A66e  ’  C12  =  ”^A12  +  A66^a6 

C13  =  °*  C14  “  ’°2b11  "  B6602,  C15  =  "^B12  +  B66^aB 

C22  3  -A22®2  "  A66a  •  C23  =  °*  C24  3  C15  ’ 

2  2  2  2 
C25  3  'B22e  "  A66a  ’  C33  =  A55  “  6  A44 

C34  3  "aA55*  C35  *  _6A44’  C44  3  "Dlla2  "  D66b2  '  A55 

C45  3  "(D12  +  W®6  ’  C55  3  "W2  "  O2202  "  A44 
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Thus,  for  a  given  a  =  mir/a,  a  =  mr/b,  q^  ,  F^"  (see  (10)  and  ( 1 2)X  and  cross 
ply  construction,  one  needs  to  solve  the  5  by  5  matrix  equation  (15)  for  the 
vector  {a}  of  amplitudes  of  the  generalized  displacements. 

FINITE-ELEMENT  FORMULATION 


As  pointed  out  in  the  previous  section,  exact  solution  to  Eq.  (8)  can 
be  obtained  only  under  special  conditions  of  geometry,  edge  conditions,  loadings, 
and  lamination.  Here  we  present  a  simple  finite-element  formulation  which  does 
not  have  any  limitations  (except  for  those  implied  in  the  formulation  of  the 
governing  equations). 

Suppose  that  the  region  K  is  subdivided  into  a  finite  number  N  of  subregions 
or  finite  elements,  Fg  (e  =  1,2,...,N).  Over  each  element  the  generalized 
displacements  (u0'v0»w»'f'x»'('y)  are  interpolated  according  to 


uQ  -  I  Vi 


r  1  s  2 
vQ  -  Z  V-  .  w  -  Z  w.*. 


(18) 


*x  s  2  *X1*1 


ib  =  Z  ib  . 


yri 


where  (a  =  1,2,3)  is  the  interpolation  function  corresponding  to  the  i-th  node 
in  the  element.  Note  that  the  in-plane  displacements,  the  transverse  dis¬ 
placement,  and  the  slope  functions  are  approximated  by  different  sets  of 
interpolation  functions.  While  this  generality  is  included  in  the  formula¬ 
tion  (to  indicate  the  fact  that  such  independent  approximations  are  possible), 

we  dispense  with  it  in  the  interest  of  simplicity  when  the  element  is  actually 

1  2  3 

programed  and  take  $.  *  =  $..  (r  =  s  =  p).  Here  r,  s,  and  p  denote  the 

number  of  degrees  of  freedom  per  each  variable.  That  is,  the  total  number  of 
degrees  of  freedom  per  element  is  2r  +  s  +  2p. 
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Substituting  Eq.  (18)  into  the  Galerkin  integrals  associated  with  the 
operator  equation  (8),  which  must  also  hold  in  each  element  Rg, 


(CL]{5>  -  dx  dy  =  0 


(19) 


and  using  integration  by  parts  once  (to  distribute  the  differentiation 
equally  between  the  terms  in  each  expression), we  obtain 


'[K11]  [K12]  [0]  [K14]  [K15]" 

fu} 

fF1* 

[K12]  [K22]  [0]  [K24]  [K25] 

fv> 

<F2} 

[0]  [0]  [K33]  [K34]  [K35] 

fw} 

= 

rf3\ 

[K14]  [K24]  [K34]  [K44]  [K45] 

f  F41 

[K15]  [K25]  [K35]  [K45]  [K55] 

f »!» v> 

fF5} 

e 

J 

e 

(20) 


where  the  {u>,  {v},  etc.  denote  the  columns  of  the  nodal  values  of  u,  v, 
respectively,  and  the  elements  (a,B  =  1,2,..., 5)  of  the  stiffness  matrix 
and  of  the  force  vector  are  given  by 

Kl]  ■  AllG?j  *A16(GU  *  +  fl66GU 

A]  ’  A12Gij  *  A16GU  *  A26G1j  +  A66G)i 

kJA  -  *  B„(Hj»  *  «*)  * 

Klj  ’  812H?)  +  816H?J  *  "atfj  *  8«>"ii 

•  Wtf  *  GJi>  *  A22G(j  *  A66G1j 

Kij  =  B16Hi j  +  B66Hi j  +  B12Hji  +  B26Hij 

*?j  -  B26(HiJ '  *  HJ?>  *  B66Hi j  *  hAl 

<33  "  A4481  j  *  A45<S?J  *  SjX(>  *  A55SU 

k34  ■  a44R?S  *  a45r>°  ,  k*  ■  a45r“  *  a55r^ 
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KU  =  D11T«  *  °16<T?j  *  ^  +  °66Ttj  *  A44T?j 
'  °26<T«  *  T£>  *  °66T«  ♦  0„TL  4  A„T° 


22  i  j  55  i j 


F?  =  |  f^]  dx  dy,  (a  =  1,2;  i  =  1,2 . r) 


-3  _ 


*1- 


qi>f  dx  dy,  (i  =  1 ,2,...  ,s) 


fa<P^  dx  dy,  (a  =  4,5;  i  =  l,2,...,r) 


where 


Gfn  = 

1J 


u^n  _ 

Hij  " 


<-^n  _ 

ij 


pCn  _ 

Rij 


ft  - 

ij 


1  1 
i  .5  *j,n 


dx  dy  ,  (i,j  =  1 ,2, . . . ,r) 


F_ 


dx  dy  ,  (i  =  1,2 . r;  j  =  1,2,. ..,t) 

*  f s  j  *n 


f  n  dx  dy  ,  ( i ,  j  =  l,2,...,s) 

p  i  J  9 n 

dx  dy  ,  (i  =  1 ,2, . . . ,s;  j  =  l,2,...,t) 

•  js  J » n 


R_ 


dx  dy  ,  (i.j  =  1 ,2, . . . ,t)  ,  (5,n  =  0,x,y) 

*  J  >T1 


(21) 


(22) 


(23) 


,xx  _  „x 
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and  =  G?j,  etc.  In  the  special  case  in  which  ^  =  <j>i  =  4>? ,  all  of  the 
matrices  in  Eq.  (23)  will  coincide. 

In  the  present  study  rectangular  elements  of  the  serendipity  family 
are  employed  with  the  same  interpolation  for  all  of  the  variables.  The 
resulting  stifness  matrices  are  20  by  20  for  4-node  element  and  40  by  40 
for  the  8-node  element.  As  pointed  out  in  a  recent  study  [6],  the  YNS 
theory  can  be  derived  from  the  corresponding  classical  thin  plate  theory 
by  treating  the  slope-displacement  relations 
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3W  _  3W 

3x  "  ”sx  ’  dy 


(24) 


as  constraints.  Indeed,  when  the  constraints  in  Eq.  (24)  are  incorporated 
into  the  classical-thin  plate  theory  by  means  of  the  penalty  method,  the 
resulting  equations  correspond  to  the  YNS  theory  with  the  correspondence, 

e x  -  i>x  and  9y  -  *  (25) 

It  is  now  well-known  that  whenever  the  penalty  method  is  used,  the  so-called 
reduced  integration  [11]  must  be  used  to  evaluate  the  matrix  coefficients 
in  Eq.  (21).  That  is,  if  the  four-node  rectangular  element  is  used,  the  1  x  1 
Gauss  rule  must  be  used  in  place  of  the  standard  2x2  Gauss  rule  to  numerically 
evaluate  the  coefficients  K. ■ .  For  more  details  on  the  effect  of  reduced 

•  J 

integration  on  the  solution  accuracy,  one  can  consult  [12,13]. 

NUMERICAL  RESULTS 


The  finite  element  based  on  the  formulation  of  previous  sections  is 
employed  in  the  static  analysis  of  rectangular  plates.  The  effect  of 
boundary  conditions,  laminations,  and  loadings  on  the  bending  deflections 
and  stresses  is  investigated.  In  the  following  analyses,  two  types  of 
materials  properties,  typical  of  advanced  fiber-reinforced  composites, 
are  employed: 

Material  I:  E^Eg  =  25,  G12/E2=  0.5,  G23/E2  =  °-2’  v12  =  °”25 
Material  II:  E^Eg  =  40,  G12/E2  =  0.6,  G23/E2  =  0.5,  v,2  =  0.25 

It  is  assumed  that  G^2  *  G23  and  vi£  =  vi3#  ^  value  of  5/6  was  used  for 

2  2 

the  shear  correction  coefficients,  »  k2  (see  Whitney  [14]).  All  of  the 
computations  were  carried  out  on  an  IBM  370/158  computer. 

First  the  effect  of  various  boundary  conditions  on  the  bending 


deflection  is  investigated.  Table  1  contains  deflections  for  six  types 
of  boundary  conditions.  BC-I  corresponds  to  a  simply  supported  case,  BC-II 
corresponds  to  the  case  in  which  two  vertical  sides  are  simply  supported 
and  the  other  two  sides  clamped,  and  BC- III  through  BC-VI  correspond  to 
various  clamped  cases.  Note  that  there  are  great  differences  between  the 
deflections  obtained  by  BC- III  and  BC-V.  The  finite  element  results  are 
in  close  agreement  with  the  exact  solutions  of  Timoshenko  [15]  and  Das  and 
Rath  [16].  Table  2  contains  bending  deflections  for  BC-I  and  BC-II  for 
isotropic  (v  =  0.3)  plate  subjected  to  uniform  temperature  distribution 
(i.e.,  q  =  0,  P.j  =  0,  -  0,  TQ  =  0,  T-j  =  1.0).  The  finite  element  results 

are  in  good  agreement  with  Timoshenko  [15]  for  all  aspect  ratios.  The 
table  also  shows  the  numerical  convergence  of  linear  and  quadratic  finite 
elements. 

Figures  la  and  lb  show  variation  of  the  bending  deflection  with  the 
aspect  ratio,  and  side-to-thickness  ratio  for  isotropic  and  orthotropic 
single-layer  plates.  Note  that  the  deflection  increases  with  decreasing 
aspect  ratio  for  thin  isotropic  plates,  whereas  for  orthotropic  plates 
(Material  I,  a 2  =  3a-j),  the  deflection  increases  slowly  with  a/b  >  1 
and  then  decreases.  That  is,  orthotropic  materials  have  a  damping  effect 
on  the  transverse  deflection.  The  effect  of  thickness  shear  strain  on  the  de¬ 
flection  is  shown  in  Fig.  lb.  for  different  aspect  ratios  of  single-layer 
orthotropic  plate.  It  is  clear  that  there  is  less  than  10%  increase  in 
deflection  for  thick  plates  (compared  to  the  thin  plate  deflection)  sub¬ 
jected  to  thermal  loading,  whereas  it  is  known  that  the  increase  in  de¬ 
flection  is  of  order  30%  in  plates  subjected  to  mechanical  loading.  Again, 
the  finite  element  results  are  in  close  agreement  with  the  closed- form 
solution  presented  herein. 
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Table  1.  Effect  of  boundary  conditions  on  the  deflection  for  isotropic 
plate  subjected  to  uniform  temperature  (q=0,  TQ=0,  v=0.3) 


1 

Mesh 

a/b 

_ 

Deflection,  w 

Boundary  Conditions 

Exact 

FEM 

1.0 

0.9578* 

0.9575 

u=w=ik,  =0  on  CD 

Q4 

1.5 

0.5824 

0.5822 

BC-I : 

2.0 

0.3702 

0.3701 

v=w=^=0  on  BC 

1.0 

0.206+ 

0.2063 

u=v=w=<|<  =0  on  CD 

Q4 

1.5 

0.036 

0.03601 

J 

BC"II;  v=w=*  =0  on  BC 

2.0 

0.0055 

0.00561 

. 

y 

1.0 

_ 

0.0137 

Q4 

1.4 

0.0224 

BC-III:  u=v=w=ii<  =0  on  CD  and  BC 

2.0 

0.0276 

1.0 

0.0138* 

0.0138 

u=v=w=ij/  =0  on  CD 

1  Q4 

1.4 

0.0226 

0.0226 

BC-IV :  y 

2.0 

0.0277 

0.0277 

u=v=w=<|j  =0  on  BC 

-  -  - - - — 1 

1.0 

0.0445 

u=v= w=^x=0  on  CD 

L4 

1.2 

— 

u=V=w='l<  =0  on  BC 

L4  | 

1.0 

— 

0.0455 

BC-IV:  u=v=w=0  on  CD  and  BC 

*  Timoshenko  [15] 


_  int  -u 

+  Das  and  Rath  [16];  w  =  w  ,  Ti  =  Tx 


Quarter  plate 


Q4  =  4x4  mesh  of  (8-node)  quadratic  elements 
L4  =  4x4  mesh  of  linear  elements 


!"/*»  01 
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Figure  2  shows  the  mechanical  response  of  two-layer  cross-ply  (0°/90°) 
square  plates  (Material  I)  subjected  to  sinusoidally  and  uniformly  distributed 
loading  (T  =  Tx  «  0).  Note  that  the  shear  deformation  effect  is  relatively 
more  pronounced  for  sinusoidal  loading  than  for  uniform  loading.  In  Table 
3,  the  deflections  due  to  thermal  loading  and  combined  loading  are  compared 
with  the  corresponding  closed -form  results.  The  table  also  contains  deflec¬ 
tions  for  single-layer  and  three-layer  (0°/90o/0°)  plates.  It  was  noted 
that  the  deflections  obtained  for  four-layer,  symmetric  cross-ply  plates 
(0°/90o/90o/0°)  are  very  close  to  those  obtained  for  single-layer  plates  (for 
the  same  total  thickness). 

The  effect  of  thickness  and  aspect  ratio  on  the  thermal  and  mechanical 
response  of  cross-ply  plates  (0°/90°,  0°/90°/0°,  0o/90°/90o/0° ,  Material  I, 

Tj  =  1.0,  a2  =  3a!,  ai  =  10"6,  Pj  =  P2  =  0)  is  shown  in  Figure  3.  The 
finite-element  solution  is  in  close  agreement  with  the  closed-form  solution 
everywhere  except  for  small  values  of  a/t  (i.e.  for  thick  plates).  The  effect 
of  thickness  shear  on  the  deflection  is  less  with  the  increasing  aspect  ratio. 
Clearly,  the  antisymmetric  cross-ply  plates  (0°/90°)  have  different  response 
characteristics  with  respect  to  the  aspect  ratio  when  compared  with  the 
symmetric  cross-ply  plates  (0°,  0°/90°/0°,  or,  0o/90°/90o/0°) . 

In  Figure  4,  the  closed-form  and  finite -element  solutions  of  four- 
layer,  cross-ply  (0°/90o/90o/0°)  plates  subjected  to  sinusoidal  temperature 
distribution  and/or  mechanical  loading  are  compared.  It  also  contains 
finite -element  solutions  for  simply  supported  platessubjected  to  uniform 
temperature  distribution,  and  clamped  platessubjected  to  parabolic  tempera¬ 
ture  distribution  along  y-axis  and  constant  along  x-axis.  No  closed-form 
solutions  are  available  for  these  two  problems.  Specifically,  the  figure 
shows  non-dimensional  deflection  versus  side-to-thickness  ratio  for  the 


(a)  Displacements 


(b)  Stresses 

Figure  2  Mechanical  bending  of  two-layer  (0°/90°) 
square  plates  (Material  I). 
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Table  3.  Effects  of  loading,  lamination,  and  thickness  on  the 
non-dimensional ized  deflections  for  simply-supported 
(BC-I)  square  plate  (Material  I,  a2  =  30^ 


a/t 

w  (sinusoidal 

temp . ) 

w  (sinusoidal  temp. 

and  loading) 

0° 

0790° 

079070° 

0° 

0790° 

079070° 

100 

1 .031 3+ 

1.6765 

1.0949 

1.0008 

2.4563 

1 .0025 

(1.0312) 

*  (1.6764) 

(1.0948) 

(1.0006) 

(2.456) 

(1 .0018) 

50 

1 .0317 

(1 .0317) 

1 .6765 

(1.6764) 

1.0963 

(1.0962) 

1.0117 

(1.0116) 

2.4625 

(2.462) 

1.01 50 

(1.0149) 

25 

1 .0334 

(1 .0333) 

1.6765 

(1.6764) 

1.1018 

(1.1017) 

1 .06068 

(1.0657) 

2.509 

(1.509) 

1.0802 

(1 .0800) 

20 

1.0346 

(1 .0345) 

1.6765 

(1.6764) 

1.1058 

(1.1057) 

1.1117 

(1.1116) 

2.5448 

(2.544) 

1.1292 

(1.1290) 

12.5 

1.0396 

(1.0395) 

1.6765 

(1.6764) 

1.1224 

(1.1223) 

1.2974 

(1.2973) 

2.7003 

(2.700) 

1.3372 

(1.3371) 

1 .0440 

1 .6765 

1.1365 

1 .4672 

2.8440 

1.5233 

i 

10 

(1.0439) 

(1.6764) 

(1.1364) 

(1.4670) 

(2.844) 

(1 .5231  ) 

6.25 

1.0602 

(1.0601) 

1.6765 

(1.6764) 

1.1870 

(1.1869) 

2.1869 

(2.1867) 

3.4667 

(3.466) 

2.2832 

(2.2829) 

c 

1.0721 

1.6765 

1.2224 

2.8332 

4.0416 

2.9424 

D 

(1.0720) 

(1.6764) 

(1.2224) 

2.8329) 

(4.041) 

(2.9421) 

+closed 

-form  solution; 

♦finite -element  solution 

w  *  10  hw/a^a2  ,  w  = 

»<r>  A  •  • 

=  ir4[4G12  + 

E1  +  (UV12?e 

1 -vi2V2i 

-  ]/12 
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following  cases (obtained  using  2x2  mesh  of  quadratic  elements): 

1.  simply-supported  square  plate  (SS)  subjected  to  sinusoidal  loading 
(SL)  (Material  I,  TQ  =  T j  -  0,  Pr  =  P2  =  0)^> 

2.  same  as  Case  1,  except  Material  II,  and  sinusoidal  temperature  (ST) 
distribution  (Tj  =  1.0)  — 

3.  same  as  Case  1,  except  sinusoidal  temperature  (ST)  distribution 
(Ti  -  102) -o_ 

4.  same  as  Case  1,  except  qQ  =  0,  and  sinusoidal  temperature  (ST) 
distribution 

5.  same  as  Case  1,  except  qQ  =  0,  and  uniform  temperature  (UT) 
distribution  — 

6.  clamped  square  plate  (CC)  subjected  to  parabolic  temperature 
distribution,  equivalent  to  the  mechanical  loading,  Pj  =  P2  =  0, 
q  =  P*,  where  — 

Ei(ai+  v2i<*2)  _ 

p*  =  Tj 

6(1 -Vj2v21 ) 

Different  nondimensional izations  are  used  for  pure  mechanical  loading  (w=  ) 

and  pure  temperature  loading  (w  =  ■^wt  ).  In  the  case  of  combined  loading, 

01 1 T  i  a2 

the  same  nondimensional ization  as  that  in  mechanical  loading  is  used.  It 
is  found  that  the  thermal  bending  (qQ  =  0)  is  virtually  independent  of  the 
mechanical  properties  (i.e.  same  for  Materials  I  and  II)  of  the  plate.  Also, 
the  thermal  bending  (w)  is  almost  independent  of  the  thickness  for  sinusoidal 
temperature  distribution.  However,  it  is  clear  from  Case  6  that  the  thickness 
effect  is  more  pronounced  on  the  deflection  of  a  clamped  plate  under  parabolic 
temperature  distribution. 

Figure  5  shows  the  effect  of  thickness  on  the  thermal  bending  (qQ  =  0, 

Pi  *  P2  5  0)  of  cross-ply  and  angle-ply  plates.  First  note  that  the  deflection 
scale  is  amplified  (compared  to  Figure  4)  in  order  to  show  the  relative 


10  20  3  0  4  0  5  0  a/t 

Figure  4  Comparison  of  closed-form  solutions  and  finite 
element  solutions  for  four-layer  cross-ply 
(0o/90o/90o/0°)  plate. 
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effect  of  thickness  shear  strain.  For  example,  if  the  deflection  for  4-layer 
(0°/90o/90o/0°)  cross-ply  (Material  II)  shown  in  Figure  5a  were  plotted 
to  the  same  scale  as  that  used  in  Figure  4,  it  would  have  overlapped  on  that 
shown  (in  Figure  4)  for  4-layer  (0°/90o/90o/0°) ,  cross-ply  (Material  I). 
Although  not  plotted,  the  deflection  vs.  side-to-thickness  ratio  plot  for  the 
four-layer  antisymmetric  cross-ply  (0°/90°/0°/90^ ,  Material  II,  a/b  =  1) 
plate  is  almost  identical  to  that  shown  for  the  4-layer,  antisymmetric  angle- 
ply  (450/-45°/45°/-450,  Material  II,  a/b  *  1)  except  for  an  additive  constant 
(i.e.  shift)  of  unity  (with  respect  to  the  nondimensionalization  used  there). 
From  Figure  5a  it  is  clear  that  the  thickness  shear  effect  is  amplified  for 
larger  aspect  ratios  (a/b).  Figure  5b  shows  the  normal  shear  stresses  for 
the  two  cases  for  which  deflections  are  plotted  in  Figure  5a.  The  thermal 
bending  (deflection  as  well  as  stresses)  for  symmetric  angle-ply  (45°/-45°/ 
-45°/45°,  Material  II,  a/b  =  1)  plate  was  found  to  be  (not  shown  here  for 
brevity)  similar  to  that  of  the  antisymmetric  angle-ply  plates  for  the  same 
material  properties,  loading,  and  boundary  conditions,  except  for  a  small 
positive  shift  in  the  deflection. 

SUMMARY  AND  CONCLUSIONS 

A  finite-element  formulation  of  equations  governing  layered  anisotropic 
composite  plates  subjected  to  mechanical  as  well  as  thermal  loading  is 
presented.  The  element  includes  the  effect  of  shear  deformation  and  involves 
five  degrees  of  freedom  (three  deflections,  and  two  slope  functions)  per 
node.  Numerical  convergence  of  linear  and  quadratic  elements  is  shown, 
and  results  are  presented  for  cross-ply  and  angle-ply  rectangular  plates 
subjected  to  sinusoidal  and  uniform  loadings;  thermal,  mechanical,  and  com¬ 
bined  loadings  are  considered.  To  check  the  finite  element  results,  a  closed- 
form  solution  is  developed  herein  for  cross-ply  rectangular  plates  subjected 


mm* 


10  20  30  40  50  a/t 

(b)  stresses  (at  Gauss  Points) 

Figure  5  Effect  of  thickness  on  the  thermal  response  of 
cross-ply  and  angle-ply  simply  supported  plates 
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to  sinusoidal  mechanical  and/or  thermal  loadings.  The  finite -element  results 
agree  very  well  with  the  closed-form  solutions.  The  maximum  error  (about  10%) 
in  deflections  and  stresses  occur  in  the  thick  plate  region  (i.e.,  for  side- 
to-thickness  ratios  smaller  than  10).  Thus,  the  finite  element  developed 
herein  is  computationally  simple  compared  to  other  plate  and  shell  elements 
used  previously  in  the  thermal  stress  analysis  of  plates.  Extension  of  the 
present  element  to  nonlinear  analysis  seems  to  be  the  next  logical  step. 

In  that  case,  the  present  element  saves  substantial  computational  costs. 
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